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Direct observation of impurity-induced magnetism in an S* = 1/2 
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Nuclear magnetic resonance and magnetization measurements were used to probe the magnetic 
features of single-crystalline Bi(Cui_a;Zna:)2P06 with < a; < 0.05 at temperatures between 2.6 
K and 300 K. The simple lineshape of the "^^P NMR signals of the pristine compound changes 
considerably for a; > and we present clear evidence for a temperature-dependent variation of the 
local magnetization close to the Zn sites. The generic nature of this observation is indicated by 
results of model calculations on appropriate spin systems of limited size employing QMC methods. 

PACS numbers: 75.10.Pq, 76.60.-k, 05.10.Cc, 75.30.Hx 



Systems with coupled spin chains also termed spin 
ladders are currently investigated as paradigms of low- 
dimensional magnets, because they provide a rich play- 
ground for studies of quantum criticality and general 
aspects of many-body physics [l|. Experiments are in- 
tended to test the predictions of model calculations or 
to explore new unanticipated physical features of these 
systems. For spin S = 1/2 antiferromagnetic (AFM) 2- 
leg ladders the ground state is a non magnetic singlet 
(^ . Upon application of a magnetic field the gap to the 
excited triplet state may be closed a.t H = He- A field in- 
duced Luttinger liquid (LL) phase was recently observed 
in this class of materials [3| and theoretical predictions 
about its behaviour could be tested quantitatively [J]. 
For H > He and at low temperatures, interladder cou- 
pling may cause the onset of magnetic order, forbidden in 
the range oi H < He by quantum fluctuations [5J . Local 
enhancement of magnetic correlations is also predicted to 
be induced by nonmagnetic entities residing on regular 
spin sites (spin holes) also for H < He and it has been 
argued that studies of the influence of such defects in spin 
arrays are both challenging and rewarding in helpingto 
reveal the generic magnetic features of such systems [61 . 

Previous studies revealed that in BiCu2P06 the S = 
1/2 Cu ions form zigzag shaped two leg ladders along the 
6-direction of the crystal lattice [7|, |8[ . The temperature 
dependence of the magnetic susceptibility xC^) indicates 
that the exchange interactions favour the formation of 
spin-singlet dimers [8[. Complications are expected be- 
cause of the zigzag geometry of the ladders [7|, which 
allows for frustration effects due to competing exchange 
interactions along the legs. Spin holes are achieved by 
replacing Cu with divalent 5' = Zn atoms. 

In previous experiments on powder samples, 
Bi(Cui_a;Zna;)2P06 was studied by magnetization [7[, 
neutron-scattering [Sj and nuclear-magnetic-resonance 
(NMR) [91] measurements. Comparisons of the data 
with results of appropriate model calculations provided 



sets of exchange interaction parameters parallel and 
perpendicular to the ladder legs. Results from zero- 
field muon-spin relaxation (/iSR) experiments were 
interpreted to indicate some kind of spin freezing at 
x-dependent, but low temperatures for samples with 
2;>0 0. 

In this Letter we report the results of measurements 
probing the static magnetization and the ■^^P NMR re- 
sponse of smgle crystals of Bi(Cui_j;Zn^)2P06. The 
NMR data for x > Q provide the first clear experimen- 
tal evidence for the formation of local impurity-induced 
magnetization and its temperature dependence in such 
systems. Results of numerical model calculations provide 
support that our observations reflect the generic nature 
of the magnetic features of such spin arrays. 




FIG. 1: (Color online) Structural subunit of BiCu2P06 show- 
ing two zigzag Cu'^^ spin chains with frustration along the leg, 
coupled via rung exchange interactions. Each ^P nucleus is 
coupled to four Cu ions via transferred hyperfine interactions 
mediated by oxygen atoms, as shown for the top P atom. 

High-quality single crystals of Bi(Cui_2,Zni:)2P06 
with X between and 0.05 were prepared as will be de- 
scribed in Ref . lid . Both phase purity and structural qual- 
ity were checked via Rietveld refinement of x-ray diffrac- 
tion data. For a; = 0, the low density of defects on the Cu 



sublatticc was confirmed by very small values of the mag- 
netic susceptibility x ^^ lo'^ temperatures. For T < 20 K 
Curie-type contributions to x(T), proportional to x, re- 
flect the partial replacement of Cu by Zn. In Fig. [1] we 
show schematically parts of the structural unit cell which 
emphasizes the zigzag array of Cu sites, the essential en- 
vironment around the P sites and the different couplings 
that were considered in our calculations. 

A first magnetic characterization of the specimens in- 
voked measurements of the magnetic susceptibility x{T) 
using a SQUID magnetometer between 4 and 300 K in ex- 
ternal fields between 0.5 and 5 T. The NMR experiments, 
covering temperatures between 2.6 and 300 K, probed the 
ensemble of '^^P nuclear spins in an external field of 7.069 
Tesla, corresponding to a resonance frequency of 121.828 
MHz. The / = 1/2 ■^^P nuclei couple to the Cu atoms 
via transferred hyperfine interaction Ahf and hence serve 
as probes of the local magnetization and its dynamics. 
Samples were aligned with their b axis parallel to the 
applied field direction. The spectra and the spin-lattice 
relaxation times Ti were obtained with standard pulse 
techniques. The typical duration of the 7r/2 rf pulses 
was 5 ^s. At all temperatures the '^^P resonance signals 
were reconstructed by sweeping the frequency, using echo 
detection and subsequent Fourier transformation. In the 
case of a; = 0, single-pulse detection was used. The relax- 
ation times Ti were obtained by analyzing the recovery 
of the magnetization after its destruction by an aperiodic 
comb of rf pulses. 

The x(T) data of our x = sample are independent 
of the external magnetic field, at least up to 5 T. They 
confirm in part previously published results [ll| , but with 
a much smaller residual x value at low temperatures (^ 
lO^'* cm^/mole Cu), reflecting the high perfection of the 
Cu sublattice in the parent-compound specimen. 

A selection of NMR signals recorded at different tem- 
peratures is shown in Fig. [2l The line widths (FWHM) 
of the individual resonances, to a fair approximation of 
Lorcntzian shape, are of the order of 10 kHz at all tem- 
peratures. This confirms the very low concentration of 
remnant paramagnetic centers, i.e., the close to complete 
singlet formation of the Cu spins. The well preserved 
periodicity of the spin arrays encourages the application 
of exact diagonalization (ED) methods in model calcula- 
tions to be discussed below. The inset of Fig. [2] displays 
the temperature dependence of the line shift K between 4 
and 300 K, reflecting the spin susceptibility XsiT). From 
the comparison of K{T) with the x(r) data we obtain 
Am = 0.50(3) T//iB. 

All our calculations to be described below are based 
on the Hamiltonian (see Fig. [T]): 
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FIG. 2: (Color online) Temperature dependence of the '^^P 
NMR spectra of BiCu2P06. Red lines are Lorentzian fits to 
the signals. The inset shows the temperature-induced shift 
of the frequency at resonance maximum (open circles). The 
dashed lines, connecting the dots in the main frame and in the 
inset, represent the result of exact diagonalization on a 10 x 2 
system, while the solid line results from a QMC calculation, 
as outlined in the text. 



With respect to the data shown in Fig. [2] we calculated 
the magnetization of a 10 x 2 spin system with the men- 
tioned ED procedure, employing periodic boundary con- 
ditions and considering NNN interactions. Due to the 
availability of experimental parameters, the numerical 
result depends merely on the choice of the exchange- 
parameter values. Employing the same parameters for 
antiferromagnetic coupling as reported in earlier work 
based on calculations on a 6 x 2 spin system [8[ , namely 
Ji^g/^B = 137.8 K, Jrung/fcB = 73.3 K and Jfr/fce = 58.4 
K, we obtain the curve shown as dashed lines both in the 
main frame and the inset of Fig. [2l The agreement with 
the experimental data is quite good at high temperatures, 
although not perfect overall. Indeed, the same calcula- 
tions on systems comprising between 6 x 2 and 10x2 spins 
revealed considerable size effects 
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singlet-triplet gap (at 7 T) was found to be overestimated 
by ^ 23% with respect to the thermodynamic limit (4.54 
meV extrapolated value). However, with 20 spins the er- 
ror reduces rapidly to 3.4%. Another type of calculation, 
using a Hamiltonian for unfrustrated but coupled ladders, 
which was promoted in Ref. [9|, results in the solid line 
also shown in the inset of Fig. [2] The parameters here 

axe Jleg/kB = Jrung/fce = 100 K, Jinterladder/^B = 10 K 

[9|, and the ^hf value quoted above. It is obvious that 
experiment and calculation do not match in the latter 
case. 

Next we consider the results of our investigation of 
samples with a; > 0. As reported before [9|, Zn impuri- 
ties on Cu sites affect the resonance signals in the form of 
considerable broadening at low temperatures. We note, 
however, additional peaks (or satellites) in the spectrum. 
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FIG. 3: (Color online) (a) ^^P NMR spectra I{v) for Bi(Cuo.99Zno.oi)2P06 at different temperatures. Note the shift of the 
satellite signal from the low to the high frequency side of the main signal upon reducing the temperature and the substantial 
broadening of the main signal below 20 K. (b) Temperature dependence of the /(i^) maxima and widths; the symbols are 
explained in the figure. In order to guarantee the visibility of the satellite signals, the main peaks were scaled to 1 (which gives 
rise to an abrupt change in color intensity at ~ 20 K). (c) The result of a QMC simulation (for details see text). 



clearly visible in Fig. |3K, where a selection of resonances 
for X = 0.01 at different temperatures is shown. The 
most intriguing feature is the relative shift of the satel- 
lite across the main resonance upon varying the temper- 
ature. For low concentrations x, the intensity ratio is 
-^satZ-^tot ~ C'ix{l — x)""^ ~ nx, with n the number of 
affected P nuclei around a Zn impurity. This ratio is 
0.04 and 0.19 for x = 0.01 and 0.05, respectively. There- 
fore we conclude that each Zn impurity changes the aver- 
age local magnetization at four neighboring P sites. The 
dominance of the coupling of each individual P nucleus 
to 4 Cu spins with S = 1/2 is consistent with the proxim- 
ity map for P in the crystal structure of BiCu2P06 and 
agrees with a claim in previous work based on symme- 
try arguments [l3|. Model calculations to be discussed 
below confirm that the satellite signal is indeed due to 
impurity- induced changes of the local magnetization. We 
note that the satellite signal is clearly observed also in 
the X = 0.05 case (not shown). The availability of single- 
crystal samples proved crucial for the direct detection 
of a well resolved impurity signal, in contrast with the 
broad and mostly featureless resonances found in previ- 
ous NMR studies on powder samples [9|, ll3[ . 

Figure [3b summarizes the temperature dependencies 
of the NMR spectrum / (i/) for x = 0.01 with respect to 
position and width. One and two Gaussians were needed 
to fit the main signal above and below 40 K, respectively. 



The onset of a second Gaussian below 40 K is most likely 
due to enhanced correlations in parts of the spin system. 
The satellite signal is again of Gaussian shape with, like 
the main line, a substantial broadening at low temper- 
atures. Quite intriguing is the temperature dependence 
of the satellite position on the frequency axis which can 
be cast into the form i^ [MHz] = a log T [K] + b, with 
a = -0.0365, and b = 122.12. 

The generic nature of these observations is indicated 
by the results of a simulation of a 2 x 100 spin array with 
a QMC stochastic series expansion (SSE) algorithm and 
2 spin holes located at 25 arbitrary positions. Since QMC 
cannot account for frustration effects, the result shown 
in Fig. [3J: is to be considered as preliminary. Instead 
of employing a time-demanding finite temperature den- 
sity matrix rcnormalization group procedure, we aimed 
at showing that the shift of the satellite peak across the 
main line is a general feature, in agreement with the on- 
set of staggered magnetization in the neighborhood of 
an unpaired spin. The exchange couplings for the un- 
frustrated (i.e. Jf,. = 0) Hamiltonian ([T]) were taken from 
Ref. 0, namely Jicg/fce = 80 K, Jmng/Jieg = 0.87. The 
average of (Si^) over 4 NNN Cu sites provides the local 
•^^P field. We note a distinct departure of the calculated 
satellite shift from the observed logarithmic temperature 
dependence at low temperatures. As preliminary den- 
sity matrix rcnormalization group calculations indicate, 
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FIG. 4: (Color online) Tj"^(T) for Bi(Cui_:,Zn:r)2P06 with 
a; = 0, 0.01 and 0.05, respectively. 



this discrepancy is most likely due to an overestimate of 
the induced local magnetization with the chosen Hamil- 
tonian. 

The results of our measurements of the spin-lattice re- 
laxation rates of all three samples are shown in Fig. 2] 
Above 20 K, T^^{T) is practically the same for all three 
samples, exhibiting an activated exp(— A/ZcbT) temper- 
ature dependence, with A/Zcb = 57.4 ± 4 K, persisting 
at least up to 200 K. The latter most likely indicates 
a negligible coupling between spin and lattice degrees 
of freedom in BiCu2P06. The gap value derived from 
T^^{T) and, according to Ref. 14 1, scaled by a factor 
Ai/rp^/A^ = 3/2, corresponds to a spin excitation gap 
Asp/k^ = 38.3 ± 2.7 K, close to values quoted in earlier 

works 0,B El- 

Considerable differences among the samples develop at 
low temperatures, however. For a; = 0, T^^(T) decreases 
further until it reaches a plateau and is practically inde- 
pendent of temperature below 5 K [16[ . For samples with 
a; > 0, T]~ (T) passes through an x-dependent minimum 
and increases again with decreasing temperature. This 
relaxation enhancement is quite prominent for x = 0.05, 
where it reaches a maximum at Tq = 3.2 K, followed by 
a considerable reduction below Tq. 

The interpretation of the observed features T^^{T) is 
not straightforward. The data in Fig. 2] indicate that 
the low-temperature NMR relaxation is completely dom- 
inated by the defect-induced uncompensated spin mo- 
ments. Thus the maximum for x = 0.05 in Fig. |4] is 
entirely due to increasing and decreasing fluctuations in 
this subsystem of spins. 

From the experimental xiT) data and the fair agree- 
ment with numerical calculations we confirm that 
BiCu2P06 is a model compound hosting a system of 
S = 1/2 spins, coupled antiferromagnetically and pe- 
riodically arranged on weakly interacting 2-leg ladders 
[7|. The previously mentioned frustration effects due to 
competing NN and NNN exchange interactions along the 
legs cannot be neglected, as already pointed out in ear- 



lier work [8[. The replacement of Cu atoms with Zn at 
the few at.% level is reflected in the appearance of satel- 
lites in the '^^P NMR spectra. We interpret them and 
their relative shifts as evidence for the presence of un- 
compensated dimers and locally enhanced magnetization 
near non-magnetic defects, respectively. The logarithmic 
temperature dependence of the satellite shift is markedly 
different from the gap-determined shape of the curve re- 
flecting the main-signal frequency. While this particu- 
lar temperature dependence is not well understood, it 
implies that the gap in the excitation spectrum of the 
defect-induced magnetization regions is much reduced, 
or even zero. 
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